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Abstract 

Duchon's method of thin plate splines defines a polyharmonic inter- 
polant to scattered data values as the minimizer of a certain integral 
functional. For transfinite interpolation, i.e. interpolation of continuous 
data prescribed on curves or hypersurfaces, Kounchev has developed the 
method of polysplines, which are piecewise polyharmonic functions of fixed 
smoothness across the given hypersurfaces and satisfy some boundary con- 
ditions. Recently, Bejancu has introduced boundary conditions of Beppo 
Levi type to construct a semi-cardinal model for polyspline interpolation 
to data on an infinite set of parallel hyperplanes. The present paper 
proves that, for periodic data on a finite set of parallel hyperplanes, the 
polyspline interpolant satisfying Beppo Levi boundary conditions is in 
fact a thin plate spline, i.e. it minimizes a Duchon type functional. 
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1 Introduction 

Thin plate spline interpolation, one of the main algorithms for multivariable 
scattered data approximation, was introduced by Duchon in [5]- [6] and has 
since generated a vast amount of research on its theory and applications (for a 
comprehensive survey, see Wendland [18]). Among all functions F : M. d — > R 
having all partial derivatives d a F of total order |a| = p in L 2 (R d ) (where 
p > d/2) and taking prescribed values at a finite (sufficiently large) number of 
fixed scattered points in M. d , the thin plate spline (or surface spline) interpolant 
to the data is defined as the unique minimizer of the integral functional 
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Duchon's theory identifies this minimizcr as a function S e C 2p - d - 1 (R d ) which 
is polyharmonic of order p, i.e. A P S (x) = with A the Laplace operator, for 
all x <E R d except the given scattered points. The computational usefulness of S 
stems from its explicit radial basis representation as a finite linear combination 
of translates of the fundamental solution of A p in R d . 

A different interpolation problem that can be formulated in the multivariable 
case is that of constructing a surface or function that matches continuous data 
prescribed on some collection of curves or hypersurfaces. In computer aided 
design, this problem is referred to as transftnite interpolation (Sabin 13J). Typ- 
ical examples include the reconstruction of 2D surfaces from level curves or 
from track data, as well as the visualization of 3D objects from scan data. A 
variational solution for transfinite interpolation along a set of curves in R 2 is 
defined by Apprato and Arcangeli pQ (see also [H Chapter X]), and its approxi- 
mate computation is obtained by smoothing finite elements. On the other hand, 
Kounchev's polyspline method treated in the monograph [9] is based on explicit 
piecewise polyharmonic functions of any number of variables. 

The present paper studies the polyspline method in a simplified setting, 
namely for periodic data given on parallel hyperplanes. Let n, p, N be positive 
integers with 2 < p < N + 1, and r := {to, . . . , t N }, where t < ti < . . . < t N 
are fixed real numbers. Let O be the closure in R™ +1 of the union of the open 
strips 

fi i :={(t,l,)eRxR B :te(t i _i,t i )}, j 6 {1,2,..., N} . (2) 

A function S : 17 — > C is called a polyspline of order p on strips determined 
by r in O if S £ (j 2p ~ 2 (0), S is polyharmonic of order p on each open strip 
f2j, j E {1, . . . , N}, and S is 27r-periodic in each of its last n variables. In 
order to determine a unique polyspline that interpolates continuous periodic 
data prescribed on the hyperplanes {tj} x R™, j G {0, 1, . . . , N}, additional 
boundary conditions are imposed on {to} x R™ and {tN} x R n . Then, provided 
that p is an even integer, Kounchev's theory Theorem 20.14] identifies the 
polyspline interpolant as minimizer of the functional 

2 

dtdy, (3) 

[to,tjv]xT™ 

subject to the boundary and interpolation conditions, where A is now the 
Laplace operator in R" +1 , and T := [— ir,ir]. 

Recently, Bejancu [3] has introduced boundary conditions of Beppo Levi 
type in the construction of semi-cardinal polyspline interpolation to data on 
the infinite set of parallel hyperplanes {j} x R", j G {0, 1, . . .}. To incorporate 
such conditions in the above setting, let i_i := — oo, tjv+i := oo, and let the 
corresponding open strips flo, &>n+i be defined as in |2]). 

Definition 1.1 Let p > 2 be a fixed integer. A function S : R n+1 — > C is called 
a Beppo Levi polyspline of order p on strips determined by r in R n+1 if the 
following conditions hold: 
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(i) SeC 2 P- 2 (R" +1 ); 



(ii) S is polyharmonic of order p on each open strip Clj, j £ {0, 1, . . . , N + 1}; 
(Hi) S is 2ir -periodic in each of its last n variables; 
(iv) S satisfies the Beppo Levi conditions 



The space of all Beppo Levi polysplines of order p on strips determined by r in 
W l+1 is denoted by S p (r, n). 

The main results of the present paper (obtained in section 3) prove that, for 
periodic data functions prescribed on {tj} x M", j £ {0, 1, ... , N}, there exists 
a unique interpolant in S p (t, n), which minimizes the Duchon type functional 



This shows that the Beppo Levi polyspline interpolant is a genuine thin plate 
spline analog for transfinite interpolation. Moreover, our variational characteri- 
zation is valid for any order p > 2, without the restriction that p > d/2 as in (fTj) 
or that p is even as in ([3]) . In the excluded case p = 1 , the results degenerate to 
the classical existence, uniqueness, and variational characterization of harmonic 
solutions to the Dirichlet problem in each strip separately. 

Note that the periodicity assumption is of importance for both theoretical 
and practical purposes. In our proofs, it will serve to simplify certain technical 
arguments that guarantee the uniqueness of the polyspline constructions. On 
the other hand, periodic polysplines are directly applicable to the problem of 
transfinite interpolation on a cylinder domain, as demonstrated, for example, in 
the case of visualization of the heart surface in medical imaging [lOj . 

The construction of periodic polysplines on strips is reduced to that of a 
family of exponential £-splines by separation of variables. To describe this 
procedure in our setting, let S G S p (r, n) and, for each t £ 1, define the 
sequence of Fourier coefficients of S with respect to its last n variables by 



where is the dot product in R™. For a fixed ( 6 Z", note that condition 

(i) above implies € £<2p-2 Further, let |£| the Euclidean norm of £ in 
IR™ and consider the ordinary differential operator 



d a S e L 2 (K x T") , V|a|=p. 



(4) 







(6) 



with null-space 



Ker£ 4 



{ 



spanle^Ve^!*,.. 
spanjl,*,...,^- 1 } , 



if e = o. 
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From the polyharmonic condition (ii), we deduce (as in [H Lemma 2.1]) 

C^(t) = 0, Wt e (tj-utj) , Vjg{0,1,...,JV + 1}. 

Also, the Beppo Levi conditions (|4]) imply (as in [3j Eq. (4)]) 

~ f Kcr — |£|) p on the interval (— oo,i ), 
S \ Ker + |£|) P on the interval (t N , oo) . 

These observations motivate the following definition. 

Definition 1.2 For a fixed £ G Z™, £/ie function s : R — > C is called a natural 
/^-spline on r «/: 

(*j s G C 2 p- 2 (K); 

f«; £ c s (t) = 0, Vt 6 (tj_i, t,-) , Vj G {0, 1, . . . , N + 1}; 

M (| " ^ (*) = 0, Vt < t , and (| + |£|) P s (t) = 0, Vt > 

TTie space of all natural C^-splines on r will be denoted by S p ^ (r). 

Remark 1.3 For £ = 0, this definition corresponds to the well-known natural 
polynomial splines of degree 2p — 1. For £ ^ 0, the use of adjoint operators 
on the two extreme intervals in condition (iii) is equivalent to the decay of the 
natural /^-spline s at ±oo. This differs from the standard natural conditions 
for Chebyshev splines (cf. Schumaker [T51 p. 396]), which employ one and the 
same natural operator on both sides. It is of interest to note that in the case 
of semi-cardinal interpolation with /^-splines treated in [3J, the splitting of the 
operator £j into its two adjoint factors on the left /right boundary pieces is a 
direct consequence of the Wiener-Hopf factorization technique. 

The above arguments show that S G S p (t, n) implies G S p ^ (t) , V£ G Z™ . 
Conversely, our results on Beppo Levi polyspline interpolation in section 3 will 
follow from the properties of the natural exponential /^-splines. The construc- 
tion and necessary analysis of the natural yCj-spline interpolation schemes are 
contained in section 2. An original contribution of this analysis is the use of 
radial basis representations to estimate the effect of the variable parameter £ on 
the size of the Lagrange functions for natural ^-spline interpolation. 

The extension of our results to the setting of data prescribed on concentric 
spheres will have to address significant technical differences. A further problem 
of interest would be to study the convergence properties of Beppo Levi polyspline 
interpolation. 

2 Natural £^-spline interpolation 

Note that the operator £j of ([6]) and the space S p ^ (r) of natural £^-splines can 
be defined not only for £ G Z", but for any £ G K". Therefore in this section we 
work with £ G K™\ {0}, the version of the results for £ = being well-known. 
We will also assume that p > 2 throughout the section. 

Our first result asserts existence and uniqueness for the problem of natural 
^-spline interpolation at the set of knots r = {to, t\, . . . , ijv}. 
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Theorem 2.1 For any set of data values {tjq, . . . , v.n} C R, there exists a 
unique natural C^-spline on t, s £ S p ^ (r), such that 

»{t J )=y j , j G {0,1,..., N}. (7) 

We will derive the proof of this theorem from the following result. 

Theorem 2.2 Let a G S p ^ (r) be any natural C%- spline on r and let ip G C p (R) 
suc/i tfwrf V (m) G L 2 (R)/Vm G {0, . . . ,p}, and ^ (tj) = 0, Vj G {0, . . . , TV}. 
T/ien the following fundamental identity' holds: 

Another direct consequence of ((SJ is the variational characterization of the nat- 
ural £^-spline s of Theorem 12.11 Although this will not be employed further, it 
may be of interest to mention it here. 

Theorem 2.3 Given the set of values {yo, . . . , t/jv} C R, let s be the unique 
natural C^-spline on r satisfying the interpolation conditions (0). If f G C p (R) 
is any other function such that f^ m ' g L 2 (R), Vm G {0, . . . ,p}, and f (tj) = yj. 
G {0, . . . , N}, then 



p 



2 



2 



dt. 



A related question at this point is whether the statements of the last two 
theorems remain true if the left natural operator (j^ — is replaced by its 
adjoint (—4i — \£,\) P ■ One way to see that the answer is positive is to go through 
the proof of Theorem 12.21 in subsection 2.1 and make the requisite changes. A 
faster way, however, is to use the next result that will also play a significant role 
in the arguments of section 3. 

Lemma 2.4 Let a G S v .^ (r) be any natural C^-spline on r and let ip G C p (R) 
such that ■0 (m) G L 2 (R), Vm G {0, . . . ,p}. Then 

/OO / 1 \ P /J \ p 

JS-'O "<"U- |S! ) W«« (9) 

where I I = — rr^ — tt is the usual binomial coefficient. 

The proof of this identity in subsection 2.2 works by showing that, after ex- 
panding the two brackets in the left-hand side, integrating term by term, and 
collecting terms of the same power of the odd power terms vanish. Since 
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replacing ■i — |£| by 4? + |£| in the left-hand side will only change the sign of 
the odd power terms, we deduce that, under the hypotheses of Lemma HOI 

£(£ + ia)V*>(| + iei)'?®* 

Therefore Theorems 12721 and [2731 also hold with f t + |£| in place of ^ - 

It can be noticed that, for a fixed £, Theorems 12.11 12.21 and 12.31 follow 
along a rather classical route in spline theory, the only novelty being due to 
the adjoint boundary operators in Definition 11.21 of the natural /^-spline. The 
classical theory, however, does not cover the dependence of the natural £^-spline 
interpolation scheme on the variable parameter £, as will be needed in section 3. 
In order to study this dependence, for each £ G Z" and j £ {0,1,..., TV}, let 
Lgj G S p ,£ (r) be the unique Lagrange function determined by the interpolation 
conditions 

L u {t 3 ) = 1 and L u {t k ) = for k G {0, 1, . . . , N} \ {j} . (10) 

Hence, if s is the natural £^-spline that interpolates the values yo, y\, ... , j/jv 
in Theorem 12. 11 we have the Lagrange formula 

N 

s(t) = J2voLui t )' V ^K. (11) 

The following theorem estimates the effect of the parameter £ on the stability 
of the Lagrange scheme (fTT|) . 

Theorem 2.5 There exists a constant Cq — Cq (p, r) > such that, for all 
j e {0, . . . , TV}, m G {0, .., 2p - 2}, and all £ e E n with |£| > \, we have: 

<a>(i+ien. v<eR - ( 12 ) 

A similar statement for m = and t £ [iojijv] was given by Kounchev [51 
Lemma 3] in the context of an £^-spline interpolation scheme with different 
boundary conditions and an even integer p. The validity of the present result 
for all orders m G {0,..,2p — 2} is of crucial importance for the analysis of 
section 3. In subsection 2.3 we will employ an original method of proof of this 
result via 'radial basis' representations of the Lagrange functions L^j. 

2.1 Proofs of Theorems [HQ E7J, and [2731 

We start by establishing the following auxiliary result. 



dt m t,J 
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Lemma 2.6 If ip £ C 1 (R) and ip' £ L 2 (R), t/ien £/iere exists > suc/i i/ia£ 

< C-0 (l + l*| 1/2 ) . VteR. 

Proof. Combining the Leibniz-Newton formula tp (t) = tp (0) + J Q '0' (it) dit with 
the Cauchy-Schwarz inequality, we obtain for t > 0: 



t x 1/2 

du 



W{t)\ < l^(0)|- 

< hmo)i+< i/2 (7 («)i 2 ^ 



|V (u)\ du 

1/2 



1/2 



The conclusion follows by making a similar estimate for t < and letting C^p 
max{|^(0)|,||Vll£2 (K) }. ■ 
Proof of Theorem 12.21 Let denote the convergent integral of ([8]) . Since 



) P a (t) = for t < t , the integration domain of can be replaced by 



the interval [t 0l oo). Using integration by parts, 



■h 



dt 



p-i 



lei ^(*) 



r/> (t)di. 



The expression in square brackets vanishes at to due to the above condition on 
fi, and it also vanishes at oo due to the exponential decay of <r at oo and the fact 
that, by Lemma [2.61 the function [4r — |£|) P tp (t) has at most an algebraic 
growth at oo. Applying similar arguments successively, 



J 6 = (-1) 



p-i 



d_ 

dt 



dt 



-Ki) cj{t)\--\i\)^{t)dt 



Since a £ S Pt £ (r), for each j £ {0, . . . , N} there exists a constant <jj such that 

p 



d ,,, 



d 
dt 



" 1 'd 
11 

p-i 



= M - + lei 



Oj-Cp-l)! (-2 I^D" e-l€l*, Vf G (M;+i), 



where £/v+i := oo. Hence 



N 



3=0 



dt 



lei uw*. 
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, N}, as well as lim^oo e"^'* 



and using the hypotheses 
«/> (t) = (by Lemma [ 



Noting that e~ l?|t 

^(*i)=0, jG{0,l,., 
it follows that = 0. ■ 

Proof of Theorem 12.11 Note that any function s £ <S P) £ (r) is uniquely 
determined as the extension on 1 of a function s G (jip-z [to,tpf] which is in 
Ker£^ on any subinterval (£j_i,tj), j G {1, ...,iV}, and which satisfies the 
following endpoint conditions: 

f^-lfD'-w lt=t = o, 

|) p S (i)| t=t „=0, 



(/(" 



e {o,i,... ! p-2}. 



(13) 

The restriction of such a function s to each subinterval (tj-i,tj), j G {1, . . . , TV}, 
is determined by 2p coefficients. Therefore imposing on s the continuity con- 
ditions of class C 2p ~ 2 at each interior knot t\, ... ,tpj_i, as well as the end- 
point conditions (fl3|) and the interpolation conditions (JT)), we obtain a system 
of (2p — 1) (N — 1) + 2 (p — 1) + iV + 1 = 2p./V linear equations for as many 
coefficients. 

To establish the nature of this system, we assume zero interpolation data: 
Vj = 0, j G {0, 1, . . . , N}, in which case the system becomes homogeneous. 
Letting s be determined by an arbitrary solution of this homogeneous system 
and a = tj) := s in ©, we obtain s (t) G Ker(^ - for t G R. Since 
s (t,-) = 0, j £ {0, 1, . . • , AT}, and N + 1 > p, we deduce s = 0. It follows that 
the above homogeneous system has only the trivial solution, which proves the 
conclusion of the theorem. ■ 

Proof of Theorem 12.31 Letting a := s and i\) := / — s in Theorem 12.21 

d_ 

dt 



Kl) S (*)(^-KI) (f(t)-s(t))dt = o, 



hence 



< 



dt 



-lei (/(*)-*(*)) 



/(*) 



eft 



S (t) 



dt. 



The inequality can become equality only if / (t) — s (t) G Ker (4z — 
tgl. Since f — s vanishes at the knots to, ... , tjy, this implies f = s. 



\z\y for 



2.2 Proof of Lemma f2T4l 

Let denote the integral on the left-hand side of ©. Expanding the two 
brackets inside this integral and collecting all the terms that have the same 
power I of |£|, we obtain 

i=o 
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where 

h,l ■■= E ( , f r ) ( I ) jT ^ (t) r/,<P-W(t)dt. 

By convention, the two binomial coefficients in the last formula are zero if / — r > 
p or r > p, respectively. The lemma will follow by showing that in (fTl)) the odd 
power terms vanish and each of the even power terms equals its correspondent 
term on the right-hand side of 

First, if I = 2k + 1, where k S {0, 1, . . . ,p — 1}, then the term of index r of 
the sum contains the same product of binomial coefficients as the term of 
index I — r. Further, for r < k, integration by parts gives the relations 



a^-^ (t) l/;(P-l+r) (t) dt 



— OO 



= (-i;r 2r / (t) v>o- r ) (t) dt. 



Indeed, the boundary terms in each integration by parts are zero due to the 
exponential decay at ±oo of a (t) and its derivatives, as well as to the algebraic 
growth at ±oo of t/> (t) and its derivatives of order at most p — 1. The latter 
growth property is a consequence of Lemma l2.6l and the hypotheses on ip. Since 
/ — 2r is odd, we deduce that the terms of indices r and I — r of the sum J^j 
cancel, hence Jf,2fc+i = 0. 

Second, if I = 2k, where k 6 {0, 1, . . . ,p}, then, for < r < k, integrating 
by parts k ~ r times yields 



OO 



cj( p - r > (t) V( p - 2fe+r ) (*) dt 



OO 

OO 

(T 

— OO 



(-l) fe_r / a (p - fe) (t) V (p " fc) (t)tft, 



where the boundary terms in each integration by parts are zero by the same 
reasons as in the previous paragraph. Combining this with a similar argument 
for k < r < 2k, it follows that J$,2k is a multiple of the integral 



CK: 



ff( p - fe ) (t) lp(P~ k ) {t) dt, 

the actual factor that multiplies the integral being the value 

H)'EH)'(/- )(*) = (* 

The last equality is a consequence of binomial expansions in the formal identity 
(1 - Xf (1 + Xf = (1 - X 2 f '. Therefore © is true. ■ 
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2.3 Proof of Theorem 1231 

For each £ £ R™\ {0}, we first introduce the integrable fundamental solution ip^ 
of the operator ©, i.e. £f<££ = S (the Dirac mass at the origin). 

Lemma 2.7 The unique integrable fundamental solution of the operator is 

M*)= ( ~| 2 t-i eH€l|t| E^iei'|^, tet, 

^ere q = %^jr, V/ G {0, . . . ,p - 1}, and 7p = (p - l)^ 1 . 

Proof. Note that an integrable function ip^ satisfies the distributional equation 
C^ip^ = Sq if and only if its Fourier transform (p£ (it) — J H e" 4 *'Vc (*) dt is 

(-If 



By Fourier inversion, 



(*) = 5- / e*"Vf (u) <fc = i / -™du. (15) 

The last integral is evaluated in Watson [T71 §6.16(1)] as 

r costu , tt 1 / 2 / iti y-i r , , m . ., 



« 2 + icr 

where K p _i is the modified Bessel function expressible in finite terms [T71 
§3.71(12)] by 



Therefore 



■-l\{p-l-l)\(2zf 



COS tu ttp If 1 1*1 1 



dli = 7 iM/oiti^P-i S Q ^' 1*1' ' 



o (W|£| 2 )" (p-l)!(2|C|) 2p - 1 ^ 

with the coefficients q from the statement. The expression of follows. ■ 
For convenience, we will work with the following 'normalization' of ip^: 

p-i 



&(t):=(-l) p lp \Z\ 2p - 1 n (t)=e-^Y,ci\Z?\t\ l , t eM.. (16) 

=o 
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For each £ £ 1"\ {0} and j £ {0, 1, ... , iV}, define the function L^- (indepen- 
dently of Theorem 12. 1[) by the 'radial basis function' (RBF) representation 

N 

L U (*) = I] W€ (* - *fc) . t G K > (17) 

fe=0 

where the set of coefficients {a^ : fc = 0, . . . , ./V} is determined uniquely in the 
next lemma from the system of interpolation conditions (jlOp . Let 

Af 5 := (tj - i fe ))^ =0 

be the matrix of this interpolation system under representation ()17[) . 

Lemma 2.8 For every f £ R™\ {0}, the matrix is positive definite, hence 
nonsingular. Accordingly, for each j £ {0, . . . , N}, the function j defined by 
|17| ) and \10\) is the Lagrange function in S p ^ (t) arising from Theorem \2.1l 

Proof. Writing the Fourier inversion formula (I15|) in the form 

-v \f\ 2p ' 1 f e itu 

we obtain, for any column vector v = (vq, . . . , vn) T £ K^ 1 , 

N N 

j=a k=o 

Further, since the functions e itfcU , for fc G {0, . . . , AT}, are linearly independent 
(their Wronskian being the multiple of a Vandermonde determinant), it follows 
that the above inequality is strict if v =/= 0. Therefore the matrix is positive 
definite and nonsingular. 

For the last part of the lemma, it is sufficient to prove that L^.j defined by 
(|17[) and (|10p is in S p ^ (r). Differentiating the above Fourier inversion formula 
shows that fa £ C 2p_2 (R), so the linear combination (fTTJ) is in C 2p ~ 2 (R). 
On the other hand, for any non-negative integer I, any knot tk £ t, and any 
parameter f, the following formula shows that the translate of an exponential 
polynomial is spanned by exponential polynomials of the same type, namely 

(t t k ) 1 e ±lf = £ a r (l,t k , |f |) tV^I*, 

r=0 

for some coefficients a r = a r (l,tk, |f|) (cf. [H]). We deduce that the form (|17p 
of L(: j also satisfies conditions (ii) and (mj of Definition 1 1.2i as required. The 
lemma is proved. ■ 
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Lemma 2.9 There exists a constant ji = /x (p, r) > swc/i t/ia£, /or any £ £ K n 
|£| > and any eigenvalue A of M^, we have A > ^Co, where Cq = ^pliy ■ 

Proof. We follow an idea from Narcowich-Ward [11] and Schaback [14]. Note 
that the diagonal entries of are all equal to p^ (0) = c . We show that Mj 
is diagonally dominant for sufficiently large |£|. Let 

p-i 

:= e^^qa 1 , u > 0, 
i=0 



so that ^ (t) = -0 (|£| \t\), Vt £ R. Since lim„_ 
depending only on p and JV such that 



i> (u) = 0, there exists 6 > 



0( u ) < .£2- 



Let ji 



6/ min |tj 
|£ — ifel > $i from which 



Vu > 6. 

Then for all £ with |£| > /i and all j ^ fc, we have 



fc=0 



Co 

2 ' 



It follows that the matrix is diagonally dominant and Gershgorin's circle 
theorem [7J p. 395] implies A > p^ (0) — p > \co-, for any eigenvalue A of M%. ■ 

Proof of Theorem [275j First we show that there exists a constant A = A (p, r) 
such that the coefficients ajk of the RBF representation (fTT|) satisfy 



|a ifc | < A, Vj,*G{0,...,JV}, V 1^1 > 1/2. 



(18) 



Note that, if a 3 := (a,jo, . . . ,ajN) T £ M Ar+1 and ej is the jth column of the 
identity matrix of order N + 1, then a.j = e_j, i.e. aj is the jth column of 



the inverse matrix M c . Hence, for |£| > and any j, fc, Lemma 12.91 implies 



|a,fc| < 1 1 a. 



J 112 



< 



m; 



where ||aj|| 2 is the Euclidean norm in 



pJV+1 



M~ c 



1 „ 2 

min | A | — Co ' 

is the induced matrix 



norm, and the minimum is taken over all eigenvalues A of M^. If i > fx, 
then JTSJ) holds with A := —. If i < /j, then for each £ £ R" with | < 
ICI — we use Cramer's rule to express each coefficient as the ratio of two 
determinants, the denominator being det Mg. The numerator is the determinant 
of the matrix obtained from by replacing its jth column with e^. Since the 
entries <p^ (tj — tk) of depend continuously on £, it follows that there exists 
a constant bo — bo (p,r) > 0, such that \aju\ < bo for all j, fc, and all £ with 

\ < 1^1 < M- Therefore (fT8|) is obtained by letting A := max |o , ^ j- 
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Since the values ip^ (t) of (fTBj) are positive and bounded above by a constant 
for all i e R and £ e R n \ {0}, from ^TTJ) and (JUJ) we deduce that dH]) holds 
for rn = 0. The corresponding estimates for the derivatives of L^j follow in the 
same way, since (|16|) implies 

p-i 

$ n) (t) = \t\ m e-MWJ2*»,l\Z\ l \t\ l , *eK, me{l,...,2p-2}, (19) 
for certain coefficients c mj ;. This completes the proof. ■ 



3 Transfinite interpolation with Beppo Levi poly- 
splines 

For any non-negative integer r, let W r (T n ) be the space of all continuous func- 
tions / : R" — > C which are 2-7r-periodic in each variable and satisfy the condition 

H/Hr ■■= E |&|(l+l£l) r <°0, 

where := J e~ % ^ ,y 'f{y)dy, for £ £ Z™, are the Fourier coefficients of /. 
In particular, W° (T n ) is the Wiener algebra of functions with absolutely con- 
vergent Fourier series. Note that < r\ < r 2 implies ||/|| r < ||/|| r , hence 
W r2 (T n ) C W ri (T n ). It is also straightforward that C k (T l ) C W (T n ) C 
C r (T"), where fc is the least integer greater than or equal to r + ^j^. 

Now recall the setting of the Introduction, in which p > 2. First, we 
state the existence and uniqueness of transfinite interpolation with Beppo Levi 
polysplines of order p on strips. 

Theorem 3.1 For any set of data functions fj £ W 2p ~ 2 (T™), j £ {0, . . . , N}, 

there exists a unique Beppo Levi polyspline S £ S p (r, n) such that 

S(t j ,y) = f j (y), j G{0,1,...,N}, y e TP. (20) 

Let B p := B p (R x T") be the space of all functions F £ C p (R n+1 ) that 
are 27r-periodic in each of the last n variables and satisfy, in usual multi-index 
notation, 

d a F £ L 2 (R x TP) , y\a\=p. 
On Bp we define the following semi-inner product and induced seminorm: 

(F,G) Bp := J ^d a FWGdtdy, VF,G£B P , (21) 

ixi» H=p 

Note that S p (r, n) C £> p for p > 2, and that H-Fljg coincides with (|5|). 
The following orthogonality result is the polyspline analog of Theorem 



13 



Theorem 3.2 Let S G S p (r, n) be any Beppo Levi polyspline as in Defini- 
tion{Tj\ IfGeB p satisfies G (tj,y) = 0, Vj G {0, 1, . . . , N}, Vy G T", then 

(S,G) Bp =0. 

This 'fundamental identity' is eventually used to establish the characterization 
of the Beppo Levi polyspline S of Thcorcm l3.1l as minimizer of the Duchon type 
seminorm ||-|| g subject to the transfinite interpolation conditions (|20[) . 

Theorem 3.3 Given the set of data functions {f 0> fi, ■ ■ ■ , f N } C W 2p - 2 (T™) ; 
let S G S p (r x T n ) be the Beppo Levi polyspline that satisfies conditions If 
F G Bp is any other function satisfying F (tj,y) = fj (y), Vj G {0, . . . , N}, Vy G 
T™, then 

\\S\\ Bp <\\F\\ Bp - 

The proofs below will employ standard notation. In particular, Z + : = 
{0, 1, . . .} and, for a multi-index f3 = {(3\, . . . , (3 n ) G Z™ , we let \(3\ := Pi + 

... + /?„,/?! := Pi\ . ../?„!, and dg := ^ . . . jg,. 
3.1 Proof of Theorem [37T1 

To construct the natural polyspline 5 with the properties stated in Theorem l3.1i 
we start from the absolutely convergent Fourier series representations 

fj (v) = E f^' v) > 3 e {0, 1, • • . ,N}, y e TP, 

where fj£, for (eZ™, are the Fourier coefficients of fj G W 2p ~ 2 (T™). For each 
£ gZ™, by Theorem 12.11 and its classical version for £ = 0, let G S p ^ (t) be 
the unique natural £^-spline on r, such that 

i G {0,1,..., TV}. 
By the Lagrange formula (fTTj) . we have 

AT 

We will prove that the function S defined by 

S (t, y) := E % (0 e ' ;<C ' a> . y G T™, tel, (23) 

possesses all the properties required by Theorem 13. II 



14 



First, we show that S e C 2p - 2 (K x TT) and S satisfies flU}. Let m G Z+, 
/3 € Z™, and a = (m, /?), such that |a| =m + |/3| < 2p — 2. For i £ [io,ijv], we 
substitute (j2"2")) into (|2"3")) and use the estimate (fT2"|) to obtain formally 



N 



\d a S(t,y)\< EE&« 



N 



fez™ j=o 

TV 

<2CoEll/ill|a 



Vy € TP. 



For t ^ [to,tjv], we use the same estimates for all terms corresponding to £ € 
Z n \ {0} in the first line of the above display, while the modulus of the term 
corresponding to £ = is bounded above by 



N 



J2\ho \d?dPL 0d (t)\ 



This expression vanishes if \(3\ > 1, or if f3 = and m > p, since Lo.j is a 
polynomial of degree at most p — 1 outside the interval [to, tjy]. In the remaining 
case (3 = and < m < p — 1 , we use the polynomial growth estimate 



d m 

— L 4 (t) 
dt m J v ' 



(l + ltr 1 -™), V*g-[t 0) tjv],Vje{0,...,iV}, 



with the constant Ci depending only on p and r. Altogether, these estimates 
show that the series ([25)1 is absolutely and uniformly convergent on compact 
sets inlxP and can be differentiated termwise up to the total order 2p — 2. 
We also obtain (|20|) from the Lagrange conditions (|10p satisfied by L^j. 

Next, to show that S is polyharmonic of order p on each of the strips fij, 
j £ {0, . . . ,N + 1}, we use the following result. 

Lemma 3.4 Suppose that a series of functions is absolutely and uniformly con- 
vergent on compact sets in the open domain fl and each term of the series is 
polyharmonic of order p on f2. Then the sum of the series is also polyharmonic 
of order p on fi. 

This is well-known in the harmonic case p = 1. For p > 2, it was proved by 
Nicolesco [12] p. 23] based on integral mean representations of polyharmonic 
functions. In our case, each term of (12311 satisfies 




St (fie 



S s (t) = 0, 



15 



for (t, y) 6 fij, j € {0, . . . , N + 1}, since S$ G S p ^ (r). Therefore S is piecewise 
polyharmonic as required. 

It remains to prove that S satisfies the Beppo Levi condition Q. Let m £ 
Z+, /3 S Z™, and a = (m, /3), such that |a| = m + = _p. It is sufficient to 
prove that the termwise partial derivative d a of series ([25)1 is a Cauchy series 
in £ 2 [to, tjv]) X T n ), which implies convergence to its sum d a S(t,y) in 
L 2 ((R\ [i ,tAr]) x T"). Note that the partial derivative d a of order p of the 
term corresponding to £ = in vanishes due to the natural conditions 
satisfied by outside the interval [io,ijv]- If £ £ Z"\{0} and < j < N, 
relations (jTHJ) and (fTT)) give the following representation of L^.j for £ < to- 

N p-l 



fc=0 (=0 

Hence, for m = and |/3| = p, Minkowski's inequality and (fT5)) imply 

1/2 



oo -/T^ 
N 



d^L u (t)e^y 



dydt 



p-l 

<ierEME c ^ 

fe=0 i=0 
JV p-l 

<^rEE c *K 

fe=0 (=0 



1/2 



(t fc -tf e 2 iei(*-**)(ft 



1/2 



Since, by changing the variable, the last integral becomes Euler's integral 

~M e -m\u du = ( 2 0! 



21+1 ■ 



we obtain 



in 



oo JT n 



Op.^ilu"-" dydt) <C 2 |Cr ? 



1/2 



for a constant C2 depending only on p and r. For m > 1 and m 
similar argument based on the RBF derivative (fH)|) gives 

2 \ 1/2 

ar^L^We*^) dydt) <C 2 \trK 



p, a 



00 JT 



after increasing C2 if necessary. Therefore, if a € Z™ +1 is any multi-index with 
\a\ = p and is any finite subset of Z™\ {0}, another application of Minkowski's 
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inequality gives 




3=0 




< 



j=0 



d a L u (t)e 



dydt 



1/2 



^E^El^ ier* 



J=0 

Since fj G W-" 2 ^ 2 (T n ) and p - ± < 2p - 2 for 2 < p, we deduce that any 
termwise partial derivative of total order p of the series (|2"5)l is a Cauchy series 
in L 2 ((— oo, to) x T"). By similar arguments, this conclusion also holds for the 
interval (tjsr, oo) in place of (-co, to), which completes the proof that S satisfies 
the Beppo-Levi condition |4]). 

The uniqueness of a natural polyspline S satisfying the conditions of The- 
orem follows directly from the uniqueness of natural ,Q-splines established 
in Theorem 12. 11 as in the proof of uniqueness of [21 Theorem 15]. ■ 

3.2 Proofs of Theorems [372] and [3731 



The proof of Theorem l3 . 21 will employ Parseval type representations of the semi- 
inner product (|2~lj) with respect to the y variable. For any F G B p and t G R, 
define the Fourier coefficients of F with respect to its last n variables by 



e- i ^F(t,y)dy, £ G II' 



(24) 



Lemma 3.5 Let F,G G B p . If m G Z+, /3 G Z™ , and a = (m,(3) satisfy 
\a\ = m+ \(3 1 = p, t/ien t/ie following identities hold with absolutely convergent 
integrals and series: 



\d a F(t,y)\ 2 dy 



E^ 



2 a 



A' 



ff*F(t,y) d°G(t,y)dy= £ (*) ^ G ? (*)> 



// 

R T™ 



R T" 



E^ 

Cez™ 



a Q ^(t,y) d<*G(t,y)dydt 



dt 1 



;F 6 (t) 



dt, 



- — Ft (t) —G ( (t)dt. 
dt m 5 v ; dt m 5 v ; 



(25) 
(26) 

(27) 
(28) 
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Proof. Differentiating with respect to t in (j!M|) , then integrating by parts with 
respect to each component of y and using the periodicity of F in y and the 
hypothesis F e C p (W l+1 ), we obtain 



where £ ^ if ^ 0. Since <9^<9 t m F (i, •) e C(T") C L 2 (T n ), the sequence 
of Fourier coefficients of df^d^F (t, ■) on the left-hand side of the above display 
is square summable and the Parseval relations (|25[) and (J^HJ) hold. Further, 
<9^d t m .F = d a F ei 2 (Rx T") shows that formula ([27]) is obtained by integrating 
(|23|) and invoking Fubini's theorem in the right-hand side. Similarly, (|2"5|) follows 
by integrating (|2l))) . the use of Fubini's theorem being allowed due to the estimate 



1— ^€ (*) -^ G H (*) 

di m iy J dt m 4 



1/2 



1/2 



< 



E^ 



2 a 



(*) 



J J\d a F(t,y)\ 2 dydt\ (J J \d a G(t,y)\ 2 dydt 



1/2 



< 00, 



VR T" 



\K T" 



which is a consequence of two Schwarz inequalities and formula (|27p . ■ 

Proof of Theorem 13.21 Due to the absolute convergence of the integrals and 
the series of formula (1281) with F := S, we have 



(S, G) B - 



J J E ^Sfof) d a G(t,y)dydt 



T „ \a\=p 



t-^t t-^t n\ \ dt m 5 V ' d.t m 5 W ' 



£GZ™ |a|=p 



where a = (m, 0), as in Lemma [3751 Next, replacing J2\ a \=p °y ^ ne double sum 
12m=oJ2\p\=p-m an d using al = m\0\ and the multinomial identity 



E 

|/3 1 — p— rn 



we obtain 



where 



(5, G) B 



E 

m=0 



P 

m 



2(p-m) 



■T— Sf (*) (*) 

dt m 4 w di m 4 v ; 
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We now claim that the hypotheses of Lemma 12.41 and Theorem 12.21 are verified 
for a := and ip := G^. Indeed, we have seen in the Introduction that 
S £ S p (r, n) implies S% £ S p ^ (t), V( £ Z". On the other hand, G £ B p implies 
5{ £ CP (R) and, by fTQ, d[ m) £ L 2 (R), Vm £ {0, . . . ,p}, V£ £ Z"\ {0}. The 

remaining condition G^(tj) = 0, Vj £ {0, ...,iV} follows from G(tj,y) = 0, 
Vj G {0, . . . , TV}, Vy G T™. Therefore Lemma 00 and Theorem O show that 

..U- 1 ^ 1 ) d ^(jr^) d ^ dt=0 > v ^ z "' 

which completes the proof. ■ 

Proof of Theorem 13.31 This follows in a standard way now by setting G := 
F - S in Theorem O to obtain (S, F - S) Bp = 0, hence < \\F - S\\ 2 B = 

\\F\\l - \\Sf B . Equality is only possible if \\F - 5|| Bp = 0, i.e. if F-S is in the 
null space of the seminorm (f2"Tj) . Since p > 2, this amounts to d a (F — S) = 
for all multi-indices a £ Z™ +1 with \a\ — p. Taking a — (0, (3) with (3 £ Z™ , it 
follows that the Fourier coefficients of F — S with respect to the last n variables 
satisfy (F^S)^ (t) = for all £ G Z"\ {0}, t £ R. Further, taking a = (p, 0), we 

find that the remaining Fourier coefficient for £ = satisfies -^(F — S) (t) = 0, 
Vt G R. Hence (F — S) (t, y) is simply a polynomial of degree at most p — 1 in 
t. Since TV > p - 1 and [F - S) {t h y) = for all j £ {0, 1, . . . , N} and y £ T", 
we obtain F - S = 0, Q.E.D. ■ 

Acknowledgements. Thanks are due to Ognyan Kounchev (Bulgarian Acade- 
my of Sciences) and Michael Johnson (Kuwait University) for stimulating dis- 
cussions on the subject of this paper. 

References 

[1] D. Apprato, R. Arcangeli, Ajustement spline le long d'un ensemble de 
courbes, Math. Model. Num. Anal. 25 (1991) 193-212. 

[2] R. Arcangeli, M.C. Lopez de Silanes, J.J. Torrens, Multidimensional Min- 
imizing Splines. Theory and Applications, Kluwer Academic Publishers, 
Boston, 2004. 

[3] A. Bejancu, Semi-cardinal polyspline interpolation with Beppo Levi bound- 
ary conditions, J. Approx. Theory 155 (2008) 52-73. 

[4] A. Bejancu, O.I. Kounchev, H. Render, Cardinal interpolation with peri- 
odic polysplines on strips, Calcolo 44 (2007) 203-217. 

[5] J. Duchon, Interpolation des fonctions de deux variables suivant le principe 
de la flexion des plaques minces, RAIRO Anal. Numer. 10 (1976) 5-12. 



19 



[6] J. Duchon, Splines minimizing rotation- invariant seminorms in Sobolcv 
spaces, in: W. Schcmpp, K. Zcllcr (Eds.), Constructive Theory of Func- 
tions of Several Variables, in: LNM, vol. 571, Springer, Berlin, 1977, pp. 
85-100. 

[7] G.H. Golub, C.F. Van Loan, Matrix Computations, third ed., The Johns 
Hopkins University Press, Baltimore, 1996. 

[8] O.I. Kounchev, Splines constructed by pieces of polyharmonic functions, 
in: P.-J. Laurent et al. (Eds.), Wavelets, Images, and Surface Fitting, A.K. 
Peters, Wellesley MA, 1994, pp. 319-326. 

[9] O.I. Kounchev, Multivariate Polysplines. Applications to Numerical and 
Wavelet Analysis, Academic Press, London, 2001. 

[10] O.I. Kounchev, M.J. Wilson, Application of PDE methods to visualization 
of heart data, in: M.J. Wilson, R.R. Martin (Eds.), Mathematics of Sur- 
faces, Proc. of the 10th IMA Intcrnat. Conf., in: LNCS, vol. 2768, Springer, 
Berlin, 2003, pp. 377-391. 

[11] F.J. Narcowich, J.D. Ward, Norms of inverses and condition numbers for 
matrices associated with scattered data, J. Approx. Theory 64 (1991) 69- 
94. 

[12] M. Nicolesco, Les Fonctions Polyharmoniques, Hermann, Paris, 1936. 

[13] M.A. Sabin, Transfinite surface interpolation, in: G. Mullincux (Ed.), The 
Mathematics of Surfaces VI, in: IMA Conf. Ser., New Ser. 58, Oxford 
University Press, 1996, pp. 517-534. 

[14] R. Schaback, Error estimates and condition numbers for radial basis func- 
tion interpolation, Adv. Comput. Math. 3 (1995) 251-264. 

[15] L.L. Schumakcr, Spline Functions: Basic Theory, third ed., Cambridge 
University Press, 2007. 

[16] M. Unser, T. Blu, Cardinal exponential splines: Part I — Theory and filter- 
ing algorithms, IEEE Trans. Signal Process. 53 (2005) 1425-1438. 

[17] G.N. Watson, A Treatise on the Theory of Bessel Functions, second ed., 
Cambridge University Press, 1995. 

[18] H. Wendland, Scattered Data Approximation, Cambridge University Press, 
2005. 



20 



